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ABSTRACT 

We report results from a deep search for redshifted Hi21cm absorption from eight 
damped Lyman-or absorbers (DLAs) detected in our earlier optical survey for DLAs to- 
wards radio-loud quasars. Hi21cm absorption was detected from the z - 2.192 DLA towards 
TXS2039+187, only the sixth case of such a detection at z > 2, while upper limits on the 
Hi 2 lcm optical depth were obtained in six other DLAs at z > 2. Our detection of Hi 2 lcm ab- 
sorption in the eighth system, at z = 2.289 towards TXS031 1+430, has been reported earlier. 
We also present high spatial resolution images of the background quasars at frequencies close 
to the redshifted Hi 2 lcm line frequency, allowing us to estimate the covering factor of each 
DLA, and then determine its spin temperature T s . For three non-detections of Hi 2 1 cm absorp- 
tion, we obtain strong lower limits on the spin temperature, T s > 790 K, similar to the bulk 
of the high-z DLA population; three other DLAs yield weak lower limits, T s > 140 - 400 K. 
However, for the two DLAs with detections of Hi 2 lcm absorption, the derived spin tempera- 
tures are both low T s = (160 + 35)x(//0.35)Kforthez = 2.192 DLA towards TXS2039+187 
and T s = (72 ± 18) x (//0.52) K for the z = 2.289 DLA towards TXS031 1+430. These are 
the first two DLAs at z > 1 with T s values comparable to those obtained in local spiral galax- 
ies. Based on the observed correlations between spin temperature and metallicity and velocity 
spread and metallicity in DLAs, we suggest that the hosts of the two absorbers are likely to 
be massive, high-metallicity galaxies. 

Key words: ISM: general - ISM: evolution - ISM: individual - radio lines: ISM 



1 INTRODUCTION 

Selected through quasar absorption surveys, and with neutral hy- 
drogen (Hi) column densities > 2 x 10 20 cm 4 , damped Lyman- 
a absorbers (DLAs) are the high-redshift analogues of gas-rich 
galaxies in the local Universe. The nature of high-z DLAs and 
their re dshift evolution ha ve long been subjects of much inter- 
est (e.g. IWolfe et alj[2005l) . Unfortunately, our information about 
these systems is mostly based on absorption studies, as it has 
proved very difficult to detect the absorber hosts in the presence 
of the bright background quasar (although see iFumagalli et alJ 
20ld: IPeroux et afluOld : iNoterdaeme et alj|2012t iKrogager et al.l 



20121) . Absorption studies only sample a narrow pencil beam 
through the intervening galaxy and hence often lend themselves 
to ambiguous interpretation. Further, the strong low-ionization 
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metal lines commonly detected in DLAs rarely provide direct 
information about local physical conditions (e.g. density, tem- 
perature, pressure) in the interstellar medium (ISM) of the ab- 
sorbers. Such information is only availabl e for a few DLAs, sys- 
tems with simple ve locity structure (e.g. IPettini & BowenlbOQll ; 



Carswell et al. 2012), molecular hydrogen (H 2 ) absorption (e.g. 
Ledoux et alj|2003l: ISrianand et alJ 2005t). neutral c arbon absorp- 
tion (e.g.lJorgenson et al.l2009l ; ICarswell et al.l201 lh . Sin* absorp- 



tion dHowk et alJ 20051) . etc. As a result, despite years of detailed 
studies, our understanding of the host galaxies of high-z DLAs re- 
mains quite limited. 

For DLAs occulting radio-loud quasars, a comparison be- 
tween the H i column density measured from the Lyman-o- absorp- 
tion profile and the H 1 21 cm equivalent width ca n be used to obtain 
the sp in temperature (T s ) of the neutral gas (e.g. lKanekar & Briggs 
2004). This provides information about the distribution of Hi in 
different phases of the neutral ISM, and is thus one of the few 
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direct tracers of local conditions in DLAs. Of course, such stud- 
ies require the assumption that the Hi column density towards 
the optical quasar is the same as that towards the radio source. 
While the latter is usually far more extended than the optical 
quasar, especially at low frequencies, < 1 GHz, very long base- 
line interferometry (VLBI) studies at frequencies similar to the red- 
shifted Hi 21cm line frequency allow one to estimate the frac- 
tion of radio flux density arising from the compact radio core, 
and thus the low-frequency covering factor / of the DLA (e.g. 
IWolfe et al.ll976l ; lBriggs & Wolfdl983l ; lKanekar et al.l2009al) . Es- 
timating the DLA spin temperature thus requires three measure- 
ments: (1) the Hi column density from the Lyman-or line, (2) the 
Hi 21 cm optical depth, and (3) the covering factor / from low- 
frequency VLBI studies (which also requires unambiguous iden- 
tification of the radio quasar core). Note that this is still based 
on the assumption that the Hi column density measured towards 
the optical quasar is the same as towards the radio quasar core 
(i.e. ove r a spatial extent of a few hundred parsec at the absorber 
redshift iKanekar et al]|2009ah . Comparisons between Hi column 
densities measured from the Lyman-or profile and from Hi 21 cm 
emission studies along similar sightlines in the Galaxy have found 
very si milar values, with a mean ratio of unity and a dispersion of 
» 10% iWakker et al.l201ll) . A similar comparison in the z ~ 0.009 
DLA towards SBS 1549+593 also yielded good agreement be- 
twe en the H i columns estimated from the two very different meth- 
ods jBowen et al.ll200ll ; IChengalur & Kanekaill2002h . 

Unfortunately, the number of high-z absorbers with spin tem- 
perature estimates is still quite small. Wh ile more than a thou- 
sand DLAs are now known at z > 2 (e.g. IProchaska et alj|2005t 
IProchaska & Wolfe! 120091 iNoterdaeme et all 120091) . there were, 
prior to this survey, only fou r DLAs with detections of H 1 2 1 cm ab- 
sorption at these redshifts (IWolfe et al.lll985l IKanekar et"al]|2006l. 
120071 : 1 Sriana nd et al. besides about twenty non-detections 

of Hi 21cm abs orption giving strong lower limits on the spin 
temperature (e.g. Carilli et al. 1996; Kane kar"& Chengaluij|2003l ; 



IKanekar et ai]|2009bl ; ISrianand et alj|2012h . The main reason for 
the lack of T s estimates in high-z DLAs is that very few DLAs are 
known towards bright radio-loud quasars. To address this issue, we 
have carried out an optical survey of 45 quasars selected on the 
basis of their low-frequency flux den sity; this resulted in the de- 
tection of eight new DLAs at z > 2 fellison et ai]|2008h . Our first 
detection of Hi 2 1 cm absorption from this sample was reported in 
lYork et all J2007h ; here, we present results from a search for red- 
shifted Hi 21 cm absorption in the full sample, which has yielded 
the first two low spin temperature DLAs at high redshifts. 



2 THE OBSERVATIONS 

Our search for redshifted H i 21 cm absorption in the eight DLAs 
of our sample (with two sightlines, towards TXS0229+230 and 
TXS2039+187, containing two DLAs; see Table [TJ was carried 
out with the Green Bank Telescope (GBT; six DLAs), the Gi- 
ant Metrewave Radio Telescope (GMRT; 1 DLA) and the West- 
erbork Synthesis Radio Telescope (WSRT; 1 DLA). The GBT ob- 
servations (proposals AGBT06B-042, AGBT07B-008, AGBT08A- 
076 and AGBT11B-221) were carried out between 2007 and 
2011, using the PF1-450 MHz and PF1-342 MHz receivers, 2 lin- 
ear polarizations and the Spectral Processor. A 2.5 MHz band, 
sub-divided into 1024 channels, was used for all sources except 
TXS0229+230, whose two DLAs were observed simultaneously, 
using two 2.5 MHz bands, each sub-divided into 512 channels. 



This yielded a velocity resolution of ~ 3 - 8 kms~', after Hanning 
smoothing and resampling. Standard position-switching, with On 
and Off times of 5 minutes each, was used for calibration, with sys- 
tem temperatures measured using a blinking noise diode. The two 
DLAs (towards TXS031 1+430 and TXS2039+187) with absorp- 
tion features near the expected line frequency were re-observed on 
multiple occasions with a bandwidth of 1.25 MHz. In both cases, 
the reality of the Hi 21 cm absorption was confirmed by the de- 
tection of the doppler shift due to the Earth's motion between the 
different runs. 

The GMRT observations of TXS 1239+376 were carried out 
in February 2011; earlier runs in 2007 and 2009 were found to 
be affected by weak radio frequency interference (RFI) and will 
not be discussed further. A total of 28 antennas were available for 
the observations, which used the GMRT Software Backend and a 
bandwidth of 2.08 MHz, sub-divided into 512 channels; this gave a 
velocity resolution of 7.6 km s , after Hanning smoothing and re- 
sampling. The calibrators 3C147 and 3C286 were used to calibrate 
the flux density scale, the system gain and the system bandpass. 

Finally, the z = 2.752 DLA toward TXS2039+187 was ob- 
served with the WSRT, in September 2008, with a bandwidth of 
2.5 MHz and 2048 channels (i.e. a velocity resolution of 1.9 kms~' 
after Hanning smoothing and resampling). The flux density scale 
and system bandpass were calibrated with observations of 3C286 
and 3C48. 

We also used the Very Long Baseline Array (VLBA) to obtain 
high-resolution images of the radio continuum emission of the six 
background quasars of our sample, to estimate the absorber cov- 
ering factor / at a frequency close to that of the redshifted Hi 
21 cm line. Five quasars were observed with the VLBA 327 MHz 
receivers in 2008 and 2009 (proposal BK153), using a bandwidth 
12 MHz, sub-divided into 32 spectral channels and with 2 polariza- 
tions, 2-bit sampling and on-source times of ss 2 hours. A strong 
fringe finder (3C454.3, 3C84, 3C147, 3C286 or 3C345) was ob- 
served during each run for bandpass calibration; phase referencing 
was not used. For TXS0620+389, we reduced an archival VLBA 
1420 MHz dataset from 2002; this had a bandwidth of 8 MHz, 
sub-divided into 16 channels, with 2 polarizations and an on-source 
time of 8 minutes. 



3 DATA ANALYSIS AND RESULTS 

All GBT data were analysed using the AIPS++ single-dish pack- 
age DISH. After initial editing to excise RFI and Spectral Proces- 
sor failures, the spectra were calibrated and averaged together to 
measure the quasar flux density. For each 10-second spectrum, a 
second-order spectral baseline was then fit to line- and RFI-free 
channels and subtracted out, during the process of calibration. The 
residual spectra were then averaged together with weights deter- 
mined by the measured system temperatures. In a few cases, a first- 
or second-order polynomial was fit to, and subtracted from, this 
average spectrum to obtain the final H i 21 cm spectrum. 

The GMRT, WSRT and VLBA data were analysed in "clas- 
sic" AIPS. Data rendered unusable by various issues (e.g. dead 
antennas, correlator problems, RFI, etc) were first edited out, and 
standard calibration procedures then used to obtain the antenna- 
based complex gains. For each source, a number of channels were 
averaged to produce a "channel-0" dataset, and a series of self- 
calibration and imaging procedures used to iteratively determine 
the antenna gains and the final image, via standard procedures. 3- 
D imaging techniques were used for the GMRT data, sub-dividing 
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Figure 1. VLBA images of the compact radio structure of the six background quasars; note that two of the quasar sightlines contain two DLAs. The quasar 
name, DLA redshift and observing frequency are indicated at the top of each panel. 



the field into 37 facets to correct for the non-coplanarity of the ar- 
ray. For the GMRT and WSRT data, the continuum image was next 
subtracted out from the calibrated visibility data, after which a first- 
order polynomial was fit to the visibility spectra on each interfer- 
ometer baseline and subtracted out. The residual U-V data were 
then shifted to the barycentric frame and imaged in all channels to 
obtain the final spectral cube. A cut through this cube at the quasar 
location yielded the Hi 21cm spectrum, after subtracting out a 
second-order baseline. The calibration steps for the VLBA data in- 
cluded ionospheric corrections and fringe-fitting for the delay rates, 
as well as the use of online measurements of the antenna gains 
and system temperatures to calibrate the flux density scale. The 
smaller number of VLBA antennas and the relatively poor iono- 
spheric stability (due to which only 5-6 antennas could typically be 
retained) meant that the final VLBA images were produced with 
phase self-calibration alone. The task JMFIT was used to measure 
the source flux density from these images, fitting an elliptical gaus- 
sian model to the radio core emission. For three sources, the entire 
VLBA emission was found to be compact, with a single gaussian 
sufficient to model the emission. Extended emission was detected 
in the remaining three sources, TXS0229+230, TXS2039+187 and 
TXS031 1+430, and a 2-gaussian model was hence used here. In 
the case of the first two, the extended radio structure is adja- 
cent to the core emission; it is hence possible that the foreground 
DLAs also cover the extended emission. Conversely, in the case of 
TXS031 1+430, there is a second radio component located 1.42" 
to the southwest of the core emission, i.e. x 11.7 kpc at the DLA 
redshift; the foreground DLA is unlikely to cover this component 



unless the absorption arises in a large disk galaxy. We also note, 
in passing, that the core of TXS031 1+430 also shows evidence for 
weak extended emission. However, the poor U-V coverage implies 
that such extended features may arise close to bright continuum 
sources due to deconvolution errors. We have hence fit a single 
gaussian component to the core emission of TXS03 11+430, with 
a second gaussian fitted to the south-western radio component. The 
VLBA images of the six background quasar cores are displayed in 
Fig-E in order of increasing DLA redshift; further, a wide- field im- 
age of TXS031 1+430, showing the southwest component, is shown 
in Fig.[2] Note that, for the three sources with a 2-component Gaus- 
sian model, the more compact source component was identified 
with the radio core. 

The VLBA results are summarized in Table [T] the columns of 
this table are: (1) the quasar name, (2) the DLA redshift, (3) the red- 
shifted H 1 2 1 cm line frequency, (4) the VLBA observing frequency 
(MHz), (5) the synthesized beam (in mas x mas), (6) the root-mean- 
square (RMS) noise (in mjy beam -1 ), (7) the integrated flux den- 
sity (S vlba ? in Jy) of the Gaussian components fit to the VLBA 
image, obtained using JMFIT, (8) the total source flux density S mt 
at the VLBA frequency, obtained either from direct measurements 
or by extrapo lating from measure ments at other freq uencies [e.g. 
at 36 5 MHz iDouglas et all 1 19961) and 1400 MHz ."Condon et al] 
1998)], (9) the deconvolved angular size of the fitted components 
(in mas x mas), (10) the spatial extent of the core emission at the 
absorber redshift (in pc x pc), and (11) the covering factor / of 
the foreground DLA, obtained by taking the ratio of the core flux 
density to the total source flux density S tot . 
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Table 1. Results from VLBA low-frequency imaging of quasars behind the DLAs of our sample. 



TXS ZDLA V21cm >"VLBA 

name MHz MHz 



Beam RMS S vlba S mt Angular size 

mas X mas mjy bm~' Jy Jy mas X mas 



Spatial extent / 
pc X pc 



0620+389 2.031 468.63 1400 13.6x5.4 5.1 
2039+187 2.192 444.99 327 122x72 11 



0.64 



0.68 
0.39 



0.81 
1.92 



9 7+ 03 x 7 +L2 



129+| x 34+8 
190+g x85+« 



82+3 x 6+'° 0.79 



1081+^ x 285% 0.35" 



1048+347 2.241 438.26 327 69x58 2.7 

0311+430 2.289 431.87 327 67x47 14 



0.40 



3.12 
1.21 



0.58 
5.96 



23+3 x gi6 +5.i 



106+2 x 55+2 
100+; x 86+^ 



192+| x 72+" 



883+j^ X458+ 



0.69 
0.52° 



0229+230 2.620 392.38 327 



71 X28 



1.7 



0.18 
0.19 



0.59 



24+ 2 4 x 0.0+2" 

41+2 x ,g+5 



195+^ x0+' 95 



0.30" 



0229+230 2.683 385.67 327 



71 X28 



1.7 



0.18 
0.19 



0.59 



24+ 2 4 x 0.0+24 

41+2 x ,g+5 



194!^ 4 x0+' 94 



0.30" 



2039+187 2.752 378.58 327 



122 X 72 



11 



0.68 
0.39 



1.92 



129+^ x 34+J 2 
190+24x85+1* 



1034+^ x 272+^ 0.35" 



1239+376 3.411 322.01 327 



71 X59 



0.34 0.62 



30+3 x 12 .7 2 



226+g x 90 + _H 0.54 



Notes: " Note that the more compact component was identified as the radio core in all five cases where a 2-Gaussian model was fit to the VLBA image. The 
covering factors for the DLAs towards TXS0229+230 and TXS2039+187 assume that only the radio core is covered. If the second VLBA component is also 
covered by the foreground absorber, the covering factors would be / = 0.56 (TXS2039+187), / = 0.63 (TXS0229+230) and / = 0.73 (TXS031 1+430). In 
the case of TXS03 1 1 +430, the second VLBA component is located 1 .42" to the south-west, corresponding to a distance of a 11.7 kpc away at the DLA 
redshift. This is unlikely to be covered by the foreground galaxy, unless it is a large disk. 



Table 2. Results of the H i 2 1 cm absorption spectroscopy and estimates of the spin temperature 



TXS 


Zem 


ZDLA 


Nhi 


V21cm 


Time 


AV 


trms 


JrdV 


/ 


T s 


name 






10 20 cm- 2 


MHz 


Hrs 


km s _I 


XlO- 3 


kms ' 




K 


0620+389 


3.46 


2.031 


2.0 ± 0.5 


468.63 


0.8 (GBT) 


3.1 


6.7 


< 0.22 


0.79 


>400 


2039+187 


3.05 


2.192 


5.0 ± 1.0 


444.99 


2.2 (GBT) 


1.6 


5.7 


0.597 ± 0.053 


0.35 


160 ± 35 


1048+347 


2.52 


2.241 


3.5 ± 0.5 


438.26 


2.5 (GBT) 


3.3 


2.5 


< 0.05 


0.69 


> 2155 


0311+430 


2.87 


2.289 


2.0 ± 0.5 


431.87 


2.0 (GBT) 


1.7 


1.7 


0.818 ±0.085 


0.52 


72 + 18 


0229+230 


3.42 


2.620 


2.0 ± 0.5 


392.38 


3.3 (GBT) 


7.5 


7.8 


< 0.25 


0.30 


> 140 


0229+230 


3.42 


2.683 


5.0 ± 1.0 


385.67 


3.3 (GBT) 


7.6 


8.8 


< 0.30 


0.30 


> 270 


2039+187 


3.05 


2.752 


5.0 ± 1.0 


378.58 


12 (WSRT) 


1.9 


5.4 


< 0.12 


0.35 


> 790 


1239+376 


3.81 


3.411 


2.5 ± 0.5 


322.01 


22 (GMRT) 


7.6 


0.8 


< 0.034 


0.54 


> 2330 



Notes: The spin temperatures of the last column assume the DLA covering factors listed in the penultimate column; see Tableland the main text for 
discussion of the DLAs towards the extended sources TXS0229+230, TXS0311+430 and TXS2039+187. 



We note that it is difficult to estimate errors on the covering 
factor as the total flux density and VLBA measurements were not 
carried out simultaneously. Further, it should be emphasized that 
the deconvolved angular sizes (and hence, the estimates of the spa- 
tial extent of the core radio emission) are upper limits, as any resid- 
ual phase errors would contribute to increasing the measured angu- 
lar size. The spatial extent of the radio core at the absorber redshift 
is < 1 kpc in all cases, significantly smaller than the size of a galaxy. 

Our estimates of the covering factor assume that the absorber 
only covers the compact VLBA core and does not cover any part 
of the extended radio emission. In the case of TXS031 1+430, 
TXS2039+187 and TXS0229+230, which have extended struc- 
ture in the VLBA images (on angular scales of ~ 0.1 - 1", or 
~ 0.8 - 8 kpc at the DLA redshift), it is possible that the fore- 



ground absorber may cover some of the extended emission, imply- 
ing a higher covering factor. Our covering factor estimates for these 
DLAs should hence be treated as lower limits. 

The DLAs at z = 2.192 towards TXS2039+187 and z = 2.289 
towards TXS031 1+430 showed strong Hi 21cm absorption fea- 
tures on multiple observing runs; the latter d etection was reported, 
and discussed in detail, by |Yorketal1 ( l2007t) . and will not be de- 
scribed here. The Hi 21 cm spectrum from the z = 2.192 DLA to- 
wards TXS2039+187 is shown in Fig.[3j this has an integrated Hi 
21 cm optical depth of 0.597 + 0.053 km s , and a velocity spread 
(between nulls) of 55 kms -1 . The peak Hi 21 cm absorption red- 
shift in this system (z ~ 2.191) lies between the metal-line and Ly-a 
absorption redshifts (zm ~ 2.192 and z me t a i s ~ 2.190), although the 
errors on the latter are large due to the low resolution of the present 
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Figure 4. The six non-detections of Hi 21 cm absorption, with optical depth (10 3 X f) plotted against barycentric frequency, in MHz. The top axis in each 
panel shows velocity (in km s -1 ), relative to the DLA redshift, while the expected H i 21 cm line frequency is indicated by the dashed line in each panel. The 
panels contain spectra for [A] the z = 2.031 DLA towards TXS0620+389 (where the shaded region indicates a frequency range affected by RFI), [B] the 
z = 2.241 DLA towards TXS1048+347, [C] the z = 2.620 DLA towards TXS0229+230, [D] the z = 2.683 DLA towards TXS0229+230, [E] the z = 2.752 
DLA towards TXS2039+187, and [F] the z. = 3.4082 DLA towards TXS 1239+376. 



optical spectrum. It would be interesting to test whether there is in- 
deed a velocity offset between the H i 21 cm and metal absorption, 
through a higher-resolution optical spectrum. 

The spectra of the remaining six DLAs are shown in Fig. [4] 
in order of increasing redshift. While weak absorption features are 
visible near the expected redshifted H i 21 cm line frequency in the 
spectra of TXS0620+389 and TXS 1239+376, these are not statisti- 
cally significant (< 3cr and « 3.5cr, respectively), on smoothing to 
a resolution comparable to the width of the feature. A Kolmogorov- 
Smirnov rank-1 test finds the spectra consistent with a normal dis- 
tribution, as would be expected from random noise. However, if one 
fits a linear baseline to the spectrum of the z = 2.031 DLA towards 
TXS0620+389 (after excluding the putative feature), and subtracts 
it out, the statistical significance of the feature increases to ss 5<t. 
Unfortunately, there is very little spectral baseline at frequencies 
above 468.7 MHz and it is hence not possible to rule out the pos- 
sibility that the apparent linear baseline arises from a ripple in the 
spectrum. We have hence conservatively chosen to quote a 3cr up- 
per limit on the Hi 21 cm optical depth of the DLA, based on the 
spectrum in Fig.[4"lA1, The other four DLAs showed no discernible 
absorption at or near the redshifted H i 21 cm line frequency. 

Sriana nd et al. I d2012h have previously reported a GMRT non- 
detection of Hi 21 cm absorption from the z = 3.411 DLA to- 
wards TXS 1239+376. However, our longer integration time with 



the GMRT has allowed us to reach a deeper sensitivity. We also 
present an updated co vering factor of / = .54 from our 327 MHz 
VLB A image, while ISrianand et al] Q2012) previously estimated 
/ = 1 from a 1 .4 GHz VLB A image. The difference is likely to be 
because compact radio cores typically have inverted spectra due to 
synchrotron self-absorption, while extended radio structures have 
steep spectra and hence contribute a larger fraction of the emission 
at lower frequencies. This emphasizes the need for VLBI estimates 
of the covering factor at frequencies close to the redshifted Hi 
21 cm line frequency. Consequently, we caution that our 1.4 GHz 
VLBA estimate of the covering factor of the z = 2.03 1 DLA to- 
wards TXS0620+389 may be an over-estimate. All other DLAs 
of our sample have their covering factors measured at 327 MHz, 
within a factor of 1.5 of the redshifted H i 21 cm line frequency. 

Our final results are summarized in Table [2] whose columns 
contain: (1) the quasar name, (2) the quasar redshift, (3) the DLA 
redshift, (4) t he Hi column den sity, Nm, from the Lyman-a ab- 
sorption line ( lEllison et al]|2008h . (5) the redshifted H i 21 cm line 
frequency, V2icm in MHz, (6) the on-source time, in hours, with the 
telescope in parenthesis, (7) the velocity resolution AV, in kms" 1 , 
(8) the root-mean-square noise (RMS) at this resolution, in optical 
depth units, (9) the integrated Hi 21 cm optical depth Jt^V, or, 
for non-detections, the 3<r upper limit on j rdV, in km s _I , ( 10) the 
covering factor /, and (11) the spin temperature T s , or, for non- 
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Figure 2. Wide-field 327 MHz VLBA image of TXS031 1+430, showing 
the radio component 1.42" to the southwest of the quasar core. 
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Figure 3. A new detection of Hi 21 cm absorption, at z = 2.192 towards 
TXS2039+187, with Hi 21 cm optical depth plotted against barycentric fre- 
quency, in MHz. The top axis shows velocity (in kms~'), relative to the 
DLA redshift, while the Hi 21 cm line frequency expected from the hya 
redshift is indicated by the vertical dashed line. Note that the offset between 
the peak Hi 21 cm and hy-a absorption redshifts (« 60 kms~') is smaller 
than the error in the hy-a redshift (a 180 kms 4 ). 



detections, the 3cr lower limit on T s , in K, computed from the ex- 
pression Nhi = 1.823 x 10 18 [T,//] / TdV. For non-detections, the 
line profile was assumed to be a Gaussian with full-width-at-half- 
maximum (FWHM)=15 kms~', similar to the widths of individ- 
ual spectral components in redshifted H 1 21 cm absorbers, with the 
RMS noise computed at a velocity resolution of ~ 15 km s -1 . This 
yields a conservative upper limit to the integrated Hi 21 cm opti- 
cal depth; note that the line FWHM for cold gas (kinetic tempera- 
ture < 200 K) due to thermal broadening is < 3 kms~'. 

For six DLAs of the sample, our non-detections of Hi 21 cm 
absorption imply that we only obtain lower limits to the spin tem- 
perature. Three of these have weak lower limits on the spin tem- 
perature, T s > (140 - 400) K, while the remaining three have 
high lower limits, T s > 790 K, similar to v alues obtained so far in 
the bulk of the high-z DLA population (e.g. lKanekar & Chengalur] 
l2003l : iKanekar et ai] |2009b ) . However, the two DLAs of the sam- 
ple with detections of Hi 21cm absorption, at z = 2.192 to- 
wards TXS2039+187 and z = 2.289 towards TXS031 1+430, both 
have low spin temperatures, T s = (160 + 35) x (//0.35) K and 
T 6 = (72 ± 18) x (//0.52) K, respectively. These are the first two 
DLAs with low spin temperatures (T s < 300 K) at z > 1. Indeed, 
the spin temperature of the z = 2.289 DLA towards TXS031 1+430 
is the lowest obtained to date for any DLA, independent of redshift! 

We note that the spin temperatures for the DLAs towards 
TXS2039+187, TXS0229+230 and TXS031 1+430 assume that the 
DLAs cover only the radio core, i.e. covering factors of / = 
0.35 (TXS2039+187), / = 0.30 (TXS0229+230) and / = 0.52 
(TXS031 1+430). If the foreground DLAs cover the entire flux 
density detected in the VLBA images, the spin temperature lim- 
its would be T s > 295 K (z = 2.620 towards TXS0229+230), 
T s > 565 (z = 2.683 DLA towards TXS0229+230) and T s > 
1265 K (z = 2.752 DLA towards TXS2039+187). The spin tem- 
perature of the z = 2.192 DLA towards TXS2039+187 would be 
T s = (255 ± 55) x C/70.56) K while that of the z = 2.289 DLA 
towards TXS031 1+430 would be T s = (100 ± 25) x (J 10.11) K. 

Finally, we cannot formally rule out the possibility that both 
VLBA components detected in the image of TXS03 1 1 +430 are ac- 
tually compact steep-spectrum lobes and that the sightline towards 
the optical QSO lies in between the two lobes. If so, the spin tem- 
perature estimate for this sightline could be incorrect, as the Hi 
column density determined towards the optical QSO need not be 
the same as that towards either of the radio lobes. 



4 DISCUSSION 

The high derived spin temperatures in high-z DLAs have been a 
puzzle ev er since the first Hi 21cm absorption stu dies of DLAs 
at z > 2 JWolfe & Davislll979l ; Iwolfe et al.lll98lh . Sightlines in 
the Milky Way and local spirals typically show low spin tem- 
peratures: more than 80% of t he sightlines through the Milky 
Way and M31 in the sample of iBraun & Walterbo si dl992l) have 
T s < 300 K. Similarly, even high-latitude Milky Way sightlines 
with Hi colu mn densities above t he DLA threshold typically have 
T s < 500 K jKanekaretal]|201lh . In the case of DLAs, both high 
and lo w spin temperatures have been found in DLAs at z < 0.7 
(e.g. lie Brun et al] 1 19971; lLane et al.lll998l: IChengalur & Kanekarl 



ng 
1." 



1999; IKanekar & Chengalurl 120011 ; lEllison et al 
trast, until the present sample, every DLA at z 
500 K , with the vast majority ha v ing T s > 1000 K dWolfe & Davis 



2012b . In con- 
1 had T s > 



1 19791 ; IWolfe et al] Il98l[ 1 19851 ; ICarilli et al] [l996; Brig gs et al 
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19971: IKanekar & Chengaiuil 120031 ; iKanekar et all I2006I . 120071 : 
Srianand et a l. 2010, 2012). 



ICurran et al, l ll200l emphasized that the high T s estimates in 
high-z DLAs could arise because the low-frequency radio emis- 
sion of the quasar is extended and hence not entirely covered by 
the foreground DLA. This is because most earlier studies of high-z 
DLAs measured the ratio T s //, rather than T s itself. However, the 
low-frequency covering factors of a sample of high-z DLAs were 
recently measured using very long baseline interferometric studies 
of the background rad io quasars, and found to be 5= 0.4 in all cases 
jKanekar et afl l2009a). Further, these authors found no difference 
between the distributions of the covering factors of the low-z and 
high-z DLA samples. This implies that low covering factors are not 
the cause of the high T s values obtained in high-z DLAs. Note that 
the VLBA measurements of the covering factor of the DLAs of the 
present sample obtained / > 0.3 in all cases. 

The simplest interpretation of the high spin temperatures ob- 
tained in high-z DLAs is that the neutral hydrogen in the absorbers 
is mostly in the warm p hase (the warm neutral medium; WNM), 
with very little cold Hi {Carilli et al.ll 19961: Ichengalur & K anekar 
I2OO0I : IKanekar & Chengaluj|2003l : IKanekar et alj|2009bh . A plau- 
sible explanation for this is a paucity o f cooling avenues in the ab - 
sorbers, due to their low metallicity dKanekar & Chengalurll200lh . 
In support of this hypothesis, a predicted anti-correlation between 
DLA metallicity and spin temperature has been recently discov- 
ered, with high T s values prevalent in low- metallicity DLAs and 
low T s values in high-metallicit v systems (jKanekar et al J [2009b; 
lEllison et ai]|2012l) . Note that lCarswell et al.l d2012h obtain a high 
WNM fraction in a low-metallicity DLA, at z = 2.076 towards 
QSO2206-199. Given the above anti-correlation, high-metallicity 
DLAs would be expected to have sizeable fractions of cold Hi, 
and hence low spin temperatures, even at high redshifts. How- 
ever, although about 10 % of DLAs at z > 2 have high metallici- 
ties ([Z/H] > -0.5; e.g. lProchaska et al.ll2007h . low spin tempera- 
tures have not hitherto been obtained at high redshifts. On the other 
hand, we note that high metallicity DLAs do indeed show a higher 
probability of detection of H2 absorption (e.g. lPetitiean et al . 2006; 
iNoterdaeme et al.ll2008l) . consistent with the hypothesis that such 
DLAs have significant cold neutral medium (CNM) fractions. 

The absorbers towards TXS2039+187 and TXS031 1+430 are 
the first two DLAs at z > 1 with low spin temperatures, T s < 300 K. 
While the d etection of H 1 21 c m absorption in the latter system was 
reported by I York et al. ld2007|X no VLBI stu dies of the background 
quasar were available then. York et al. (2007]) hence used a GMRT 
602 MHz image to estimate a covering factor of « 1. Our VLBA 
observations yield the lower limit / > 0.52 on the covering factor, 
assuming that the core alone is covered by the foreground absorber, 
and confirm the low estimated spin temperature in this DLA. We 
note that the spin temperature estimate in this DLA i s lower by a 
factor of s 2 than the original estimate of lYork et al.l d2007l) . em- 
phasizing the need for VLBA imaging studies at frequencies close 
to the redshifted H 1 21 cm line frequency. 

The low T s values in the DLAs towards TXS2039+187 and 
TXS031 1+430 indicate that the two absorbers are likely to con- 
tain significant fractions of the cold neutral medium, > 50% of the 
total Hi content (as T s is the column-density weighted harmonic 
mea n of spin temperatures i n different phases along the sightline; 
e.g. IKanekar & Briggsll2004h . The anti-correlation between spin 
tempe rature and metallicity jKanekar et al.l l2009bl; lEllison et al.l 
l2012h implies that both these DLAs should have high metallicities, 
[Z/H] > -0.5; conversely, the DLAs with high spin temperatures 
would be expected to have low metallicities, [Z/H] < -1. As ex- 



pected, the absorber towards TXS031 1+430 has [Si/ H] > -0.6, 
among the highest metal licities of DLAs at z > 2 jYork et all 
l2007l : lEllison et al.l |2008), while, in the case of the system to- 
wards TXS 2039+187. lEllison et all d2008h obtain [Si/H] > -1.4. 
As noted bv lEllison et alj J2008h . both of these estimates are lower 
limits as they are based on low-resolution (« 5A) spectra and as- 
sume that the Sin-yl 1 808 absorption lines are unsaturated. It would 
be interesting to accurately determine the metallicity of the DLA 
towards TXS2039+187 from a high-resolution spectrum. 

If the two DLAs do indeed have high metallicities, the cor- 
relation j3etwerai_j^ocityspj^ad^ in high-z 
DLAs jLedoux et al] 2006; Prochaska et al. 2008), interpreted as 
arising due to an underlying mass-metallicity relation, suggests 
that the two absorbers are likely to be massive galaxies. Note that 
the velocity spread (between nulls) of the Hi 21 cm absorption in 
TXS031 1+430 is indeed quite large, 130 kms~', as would be ex- 
pected of a large galaxy; it is thus plausible that the absorber also 
covers the south-west VLBA component. If so, it might be pos- 
sible to map the Hi 21cm absorption against the VLBA struc- 
ture and directly estimate the transverse extent of th e DLA (e.g. 
IKanekar & Briggsll2004l : lKanekar & Chengaludl2005l) . Conversely, 
the Hi 21cm absorption towards TXS2039+187 has a velocity 
spread of only 55 kms~'; if this is a massive galaxy, it is likely 
to be a relatively face-on system. Follow-up imaging studies of the 
host galaxies of the two DLAs will be of much interest. 

In summary, we report a deep search for redshifted H 1 21 cm 
absorption in eight DLAs at z > 2, with the GBT, the GMRT and 
the WSRT, allied with low-frequency VLBA imaging of the quasar 
radio emission to determine the DLA covering factors. High spin 
temperature limits, T s > 790 K, were obtained in three absorbers, 
and weak lower limits, T s > 140 - 400 K, in three systems. We 
report a new detection of Hi 21 cm absorption in the z = 2.192 
DLA towards TXS2039+187, only the sixth case of Hi 21 cm ab- 
sorption in DLAs at z > 2. The detection of Hi 21 cm absorption in 
the eighth DLA, at z = 2.289 towards TXS031 1+430, was earlier 
reported by lYork et all d2007l) . but our new low-frequency cover- 
ing factor measurement allows a better estimate of the DLA spin 
temperature. These are the first two DLAs at z > 1 with low spin 
temperatures, T s < 300 K. The fraction of DLAs with low spin 
temperatures (« 7%) at z > 2 is similar to that of DLAs with high 
metallicity, [Z/H] > -0.5 (~ 8%). This supports the hypothesis that 
the high spin temperatures of most high-z DLAs arise because the 
low absorber metallicities imply fewer routes for gas cooling. 
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